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ABSTRACT 

Context. Wide-field narrowband surveys are an efficient way of searching large volumes of high-redshift space for distant galaxies. 
Aims. We describe the Wide Field Imager Lyman-Alpha Search (WFILAS) over 0.74 sq. degree for bright emission-line galaxies at z ~ 5.7. 
Methods. WFILAS uses deep images taken with the Wide Field Imager (WFI) on the ESO/MPI 2.2 m telescope in three narrowband (70 A), 
one encompassing intermediate band (220 A) and two broadband filters, B and R. We use the novel technique of an encompassing intermediate 
band filter to exclude false detections. Images taken with broadband B and R filters are used to remove low redshift galaxies from our sample. 
Results. We present a sample of seven Lya emitting galaxy candidates, two of which are spectroscopically confirmed. Compared to other 
surveys all our candidates are bright, the results of this survey complements other narrowband surveys at this redshift. Most of our candidates 
are in the regime of bright luminosities, beyond the reach of less voluminous surveys. Adding our candidates to those of another survey 
increases the derived luminosity density by ~30%. We also find potential clustering in the Chandra Deep Field South, supporting overdensities 
discovered by other surveys. Based on a F0RS2/VLT spectrum we additionally present the analysis of the second confirmed Lya emitting 
galaxy in our sample. We find that it is the brightest Lyor emitting galaxy (1 x 10"'^ ergs"' cm"^) at this redshift to date and the second 
confirmed candidate of our survey. Both objects exhibit the presence of a possible second Lyo- component redward of the line. 
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1. Introduction 

Dete ctions of both galaxies and QS Os at z ~ 6 jpan et all 

E002 ; Becker et al. 2001; Diorgovsk i'et alJl200lh indicate that 
the Universe was largely reionised at that epoch. The re- 
cent three-year WMAP results combined with other cosmo- 
logi cal surveys stiggest an epoch of reionisation around z ~ 
10 dSpergel et alJ |2003) . consistent with both QSO results 
jpan et alJl2qoir and the epoch predicted by st ructure forma- 
tion models dOnedin & Ostrikeifll997tlHaiman & Loeb 1998). 



* Based on observations made with ESO Telescopes at the La Silla 
Observatory (Programmes 67.A-0063, 68.A-0363 and 69.A-0314). 



While the UV contributions of QSOs and AGN are almo st cer- 
tainly not responsible for reionisation jBarger et alj2003t) . faint 
star forming galaxies nee d to exist in extraordinar y numbers 
if they are to be the cause jYan & Windhorstll2004 . However, 
analyses of the Hubble Ultra Deep Field failed to find sufficient 
numbers of faint galaxies to support this idea (iBunker et alJ 
l2004 iBouwens et alJEoOsI) . Therefore, it is crucial to inves- 
tigate what the contribution to the ionising UV flux is from 
young stellar populations of star forming galaxies. 

Broadly speaking, two classes of star-forming galaxy dom- 
inate high redshift surveys: Lyman Break Galaxies (LBGs) and 
Lyman-a Emitters (LAEs). LBG surveys, which now number 
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in the thousands of objects at z = 3 to 5, find clumpy source 
distributions and a two-point angular correlation function in- 
dicative of strong clustering teiavalisco & Dickinsonl '2001 ; 
Foucaud et al. 2003.: Ade l berger et al. 200 3: Ou chi et al 2004 ; 
Hildebrandtet alJ2005UAllen et alJ2005l) . LAEs also show ev- 



idence for clustering although many of the LAE surveys target 
fields surrounding known sources such as proto-clusters, radio 
galax i es and OSOs (e.s. Isteidel et alJbOOOt lumer & Fvnbc 

■ Mom 



2001 



IStiavelU et al.ll200U IVenemans et al.ll2002t lOuchi et al 
20051) . On average, LAEs number l.SxlO^deg"^ per unit 
redsh ift down to 1.5 x 10~'^ erg s~' cm~^ at z = 3.4 and 
4.5 JHu et al.l Il998h . Also, their consistently small size 
(<0.6/z"' kpc) suggests they are subgalactic clumps resid- 
ing in the wind-driven outflows of larger unseen hosts (e.g. 
iBland-Hawthorn & Nulsenll2004 . Such mechanisms provide 
a straightforward means of UV photon escape from the host 
galaxy, efficiently reionising the surrounding IGM in a way 
than ordinary LBGs can not. 

The most efficient way to find LAEs is through imaging 
surveys using a combination of broad- and narrowband fil- 
ters. The advent of wide field cameras has allowed system- 
atic imaging searches that have been carried o ut to build up 
samp l es of candidate LAEs at high redshifts (e.g.lRhoads et al.l 



l2003UAiiki et alJl2003tlH7i et aT.i2004HWang et al .1120051) . The 

availability of high throughput spectrographs on 8 to 10 m- 
class telescopes has enabled the spectroscopic confirmation of 
these galaxies. Such direct imaging searches typically cover 
10^-10"^ times the vol ume of blind long -slit spectroscopic 
searches (e.g. Table 4 in lSantos et al.l2004r Furthermore, can- 
didates from narrowband surveys always have an identifiable 
emission feature that is well separated from sky lines courtesy 
of the filter design. This is in contrast to other methods, in- 
cludin g the widely-used "dropout" technique (e.g. lSteidel et al.l 
Il99qh . 

The narrowband filter design leads to a higher candidate 
LAE selection efficiency than other techniques. The only way 
to secure the identification of the emission line is spectro- 
scopic follow-up. The most common low redshift interlop- 
ers are the emission line doublets of [On] /1/13726,3728 and 
[O m] /1/14959,5007. These can be identified by obtaining spec- 
tra with a resolution R > 1500 to separate the line pair. 
Other emission lines, such as Ho- and H/3, can be identified by 
neighbouring lines. The narrowband technique has been suc- 
cessfully applied by man y authors in order to discover galax- 
ies at redshift 5-6 (e.g. 



y authors in order to discover galax- 
lAiiki et al.ll2003t iMaier et al.ll2003t 



iRhoads et alJl2003UDawson et al.l2004tlF[u et alJ2004l) and to 
locate galaxies at redshift 6 -7 dCubv et alj2003UKodaira et alJ 
l2003t IStanwav et al.ll2004 . Likewise, we employ the narrow- 
band technique in the Wide Field Imager Lyman- Alpha Search 
(WFILAS) to find g alaxies at z ~ 5.7. In Paper I in this series 
I Westra et alJ2005h . we described a compact LAE at z = 5.721 
discovered by our survey. 

In this Paper, we describe the survey design and sample 
analysis of WFILAS. In Sect. |2lwe describe the scope of the 
survey and the observing strategy. The data reduction is de- 
scribed in Sect. |3] Section 0] outlines the candidate selection 
and Sect.|5loutlines sample properties and comparison to other 
surveys. We discuss the spectroscopic follow-up of two can- 



100 - (a) 
80 

60 - /' 

- / 
40 I 

20 ^ 
: / 
tJ^ 




4000 




7800 7900 8000 8100 8200 8300 8400 8500 
Wavelength (A) 

Fig. 1. Filter set used for the WFILAS Survey, a) The broad- 
band BR and intermediate Mgis (8 15/22 nm) passbands. b) 
The M8I5 intermediate passband (solid line) and three purpose- 
built narrowbands A^gio (810/7), A^8i7 (817/7) and Nfa4 (824/7), 
shown over the wavelength region targeted for Lyo' line detec- 
tions. The transmission curves of the narrowband filters are for 
illustrative purposes only. The OH night sky line background is 
also shown. 

didates in Sect. |6l Throughout this paper we assume a flat 
Universe with (n,j,,f2A) - (0.3,0.7) and a Hubble constant 
Ho = 7 0kms~' Mpc~'. Al l quoted magnitudes are in the AB 
system JOke & GunnI 19831) '. 

2. WFILAS Survey Design and Observations 

The sky area surveyed by the WFILAS is ~0.74sq. degree. We 
observed three fields in broadbands B, R and in an intermedi- 
ate width filter centred at 815 nm encompassing three narrow- 
band filters (Fig.0. The adoption of an additional intermediate 
width filter encompassing the multiple narrowband width fil- 
ters is a novel approach compared to previous narrowband sur- 
veys. The application of the intermediate band filter enables us 
to drastically reduce the number of spurious detections in the 
narrowband filters. The narrow width of the narrowband filters 
(FWHM = 7 nm) gives a prominent appearance to emission line 
objects. Furthermore, the three chosen fields are spread across 
the sky to enable us to average out variations in cosmic vari- 
ance. Our search has covered one of the largest co-moving vol- 
umes compared to other surveys. Table ^ compares WFILAS 
with other published surveys. 

Th e observations we re taken with the Wide Field Imager 
(WFI; iBaade et al.lfl999l) on the ESO/MPI 2.2 m telescope at 
the Cerro La Silla Observatory, Chile. The data were taken over 
65 separate nights from 2001 January 19 to 2003 December 1. 
The WFI is a mosaic of eight (4 X 2) 2k x 4k CCDs arranged 
to give a field of view of 34' x 33'. The pixels are 0'.'238 on a 
side. 

As WFILAS was planned as joint project of ESO Santiago 
and the COMBO- 17 team at MPIA Heidelberg, three fields 

' niAB = -2.5 log /,, - 48.590, where ;;?ab is the AB magnitude and 
/,, is the flux density in ergs s"' cm"-^ Hz"' 
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Survey 


Fields 


Total Area 


NaiTowband 


Filter 


Co-moving 


Narrowband Detection 






(sq. degree) 


Filters 


Width (A) 


Volume (Mpc'') 


Limit (yuJy) 


LALA iRhoads & MalhotrJ200lh 


1 


0.19 


2 


75 


0.2x10'' 


0.41 


CADIS JMaier et alJ200.3h 


4 


0.11 


8-9" 


20 


0.04x10** 


3.33 


A03 jAiiki et alj2003h 


1 


0.26 


1 


120 


0.2x10* 


0.14 


SSA22 JHu et alJ2004 


1 


0.19 


1 


120 


0.2x10** 


0.30 


WFILAS (this paper) 


3 


0.74 


3* 


70 


1.0x10* 


1.06-1.74 



Table 1. Narrowband surveys for Lya at z = 5.7 



" CADIS is based on imaging with a tunable Fabry-Perot interferometer scanning at equally spaced wavelength steps iHippelein et all2003h . 
* An additional encompassing mediumband filter was used here. 



were selected to overlap with the COMBO-17 survey, i.e. 
their extended Chandra Deep Field South (CDFS), SGP (South 
Galactic Pole) and Sll fields. The coordinates of the field cen- 
tres and the exposure times in each of the filters for each field 
are given in Table |2 All three fields are at high Galactic lati- 
tude (\b\ > 54°) and have e xtinctions less than E(B - V) = 0.022 
mag jSchlegel et al.ll998h . 

We employ standard broadband B and R filters. The inter- 
mediate band (FWHM = 22 nm) observatory filter is centred at 
815 nm. The three custom made narrowband {FWHM = 7 nm) 
filters are centred at 8 1 nm, 817 nm and 824 nm. The transmis- 
sion profiles of the filters are shown in Fig.^ The intermediate 
and narrowband filters are designed to fit in the atmospheric 
815 nm OH-airglow window, where the brightness of the sky 
background is low and hence favourable to detect Lyo- emission 
at redshift -5.7. The data taken with the intermediate band fil- 
ter confirm detections of the Lye line in one of the narrowband 
filters. The broadban d B and R d ata, which were taken from the 
COMBO-17 survey dWolf et aljEb04j) . are used to confirm the 
absence of continuum blueward of the Lya line and to avoid 
sample contamination by lower redshift emission line galaxies 
(e.g. Ha at z ~ 0.24, or [O ii] at z ~ 1 .2). 

To establish the photometric zero-point of the intermediate 
and narrowband filters two spectrophotometric standard stars 
(LTT3218 and LTT7987;.BesselI 1999J were observed. 

Between 10-50 exposures were taken for each intermediate 
and narrowband filter for each field. The exposure times varied 
between 1000 and 1800 sec per frame, with a typical exposure 
time of around 1600 sec. All frames are background-limited de- 
spite the low night sky emission in this spectral region. The 
median, first and last decile of both seeing and background are 
given in Tabled 

3. Data Reduction 

The data were processed with standard IRAF^ routines 
(MSCRED TASK) and our own specially designed scripts. 
The initial steps in the reduction process consist of removing 
the zero level offset with bias frames, normalising pixel-to- 
pixel sensitivity differences with twilight flatfield frames and 

- IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



Filter No. of Background (yuJy/n") Seeing (") 

Frames 10% 50% 90% 10% 50% 90% 



A'sio 


92 


17 


27 


36 


0.65 


0.79 


1.12 


Nsn 


75 


19 


30 


41 


0.64 


0.84 


1.16 


Ns24 


77 


17 


27 


36 


0.63 


0.80 


1.10 


Msi5 


80 


17 


22 


33 


0.65 


0.83 


1.09 



Table 3. The median, first and last decile of background and 
seeing for the WFILAS narrowband imaging for all three fields 
combined. One pixel corresponds to 0'.'238. 



removal of fringes with fringe frames. During these steps, the 
8 CCDs that make up a single WFl image are treated indepen- 
dently. These processes are described in detail below. 

Normally, the overscan region of the science frames can 
be used to remove the zero level offset. However, it was no- 
ticed that the bias frames contained significant intermediate 
scale structure (10-30 pixels). To remove this, bias frames were 
taken on every day of our observations and averaged into a bias 
frame for that day. In order to minimise the noise added to the 
data by subtracting the bias, the bias frames were smoothed 
by 5 pixels and 30 pixels in horizontal and vertical direction of 
the CCDs, respectively, and subsequently medianed. The struc- 
tures are stable over periods of several months. Therefore, it 
was possible to use bias frames from different nights without 
degrading the quality of the data. 

Typically, five twilight flatfield frames were taken in one 
night for one or more filters. The frames were medianed and 
the science data was divided by the median. Hence pixel-to- 
pixel sensitivity differences were removed. The structure in the 
individual flatfield frames was stable over a period of several 
weeks. Frames taken on different nights could thus be reused. 
Any differences between flatfield frames were due to the ap- 
pearance or disappearance of dust features, or large scale illu- 
mination differences. The differences rarely amounted to more 
than a few percent. 

The raw data in the intermediate and narrowband filters 
show fringe patterns with amplitudes of up to 10% which was 
only partially removed after the data had been flatfielded. To 
entirely remove the fringe pattern, we subtracted a fringe frame 
created from 10-30 science frames. The fringing is very stable 
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Filter 


Passband/ 


CDFS field 




Sll field 






SGP field 




r WHM 




32"" 25^34 


11' 


42"' 59';933 


00 


'' 45'" 55!024 




(nm) 


-IT 48' 49'.'75 


-or 42' 46'.' 44 


-29° 34' 55'.'05 






{a) 


(b) 


(c) 


{a) 


{b) 


(£■) 


{a) 


(.b) (c) 


Narrowband A'gio 


810/7 


48.0 


0.57 


0.79 


AAA 


0.55 


0.80 


31.5 


0.87 1.03 


Narrowband Ngn 


817/7 


41.1 


0.55 


0.79 


79.9 


0.53 


0.92 


0.0 




Narrowband A'824 


824/7 


41.0 


0.72 


0.80 


43.5 


0.81 


0.87 


42.8 


0.62 0.89 


Mediumband Mm" 


815/20 


52.7 


0.29 


0.85 


33.3 


0.38 


0.88 


18.9 


0.41 0.90 


Broadband B" 


458/97 


5.0 


0.07 


1.09 


9.4 


0.07 


0.98 


10.0 


0.14 1.22 


Broadband R" 


648/160 


15.1 


0.05 


0.75 


21.2 


0.07 


0.75 


21.5 


0.07 0.76 



Table 2. WFILAS fields, filter set exposure times and detection limits. The entries under each field heading list: (a) the total 
exposure time (ks), (b) the flux for a 2cr detection on 6 pixel diameter aperture (yuJy) and (c) the final seeing ("), in each filter. 



" Broadband B and R and part of the intermediate band Mgis taken from the COMBO- 17 survey iWolf etall2004h 



over time, so we were able to use data spanning several months. 
Certain science frames still show fringe patterns because they 
are contaminated by either moonlight or twilight. Residual dif- 
ferences in the level of the background between the different 
CCDs were removed by subtracting the median background 
level from each CCD. 

To produce the final deep images we only used images 
with a seeing of less than 5 pixels (= 1'.'2) and without signif- 
icant residual fringing. To make the combining of the images 
possible, we had to apply an astrometric correction based on 
stars from the USNO CCD Astrograph Catalogue 2 (UCAC2; 
IZacharias et al ] |2004 in the three observed fields. The frames 
have a set pixel scale of (7/238 pixel ' with North up and East 
left. The images were weighted according to their exposure 
time and combined using the IRAF "mscstack" routine reject- 
ing deviant pixels. Table|2lsummarises the depth, image quality 
and total exposure time, for each coadded frame. 

4. Sample Selection and Completeness 

4.1. Photometry and Noise Characteristics 

Initial source catalogues were created for each of the 8 narrow- 
band images. Each catalogue contains the photometry for the 
sources in all 6 filters. We used the SExtract or source detection 
software (version 2.3.2, double image mode: lBertin & Arnoutj 
Il996l) . Sources were selected when at least 5 pixels were O.Scr 
above the noise level in the narrowband image used for detec- 
tion. The photometry was measured in two apertures, 6 and 10 
pixels in diameter (= r/4 and 2''4, respectively). The 6 pixel 
aperture was used to maximise the signal-to-noise of the flux 
of the objects, while the larger 10 pixel aperture was used for 
the more accurate determination of the total flux and hence the 
star formation rate. 

Some authors have found that SExtr a ctor underes timates 
flux uncertainties jFeldmeier et al.ll2002t iLabbe etaI.l E003l 
SExtractor estimates the uncertainties using various assump- 
tions that are often not valid (e.g. perfect flatfielding, perfect 
sky subtraction). The pixel-to-pixel noise in our data is slightly 
correlated because the scatter in the counts summed in 6 pixel 
apertures is about 10% higher than what one would derive from 
the measured pixel-to-pixel RMS. 



We devised a method to correct the uncertainties given by 
SExtractor to their true values as follows. First, sources with 
flux in all filters and their /Wsis magnitude between 16 and 
23 were selected. Sources brighter than /V/gis = 16 are typi- 
cally saturated, while those fainter than M^s = 23 are incom- 
plete (see Sect. 14.31 for a further discussion of incomplete- 
ness). The Mjiis-N colour (where is any of narrowband fil- 
ters /Vgio, /Vgn, or /V824) is the same for any flat continuum 
source. Therefore, the spread in the M^ns-N colour will be 
the same as the true flux uncertainty from the two contribut- 
ing filters. Next, the sources were binned into 200-source bins 
based on their /Wsis magnitude. In Fig.|3]we plot the M^i^-N 
colour versus the /Wgis magnitude of one of our Sll catalogues. 
Mean values for the Mms-N colour, /Wsis, N magnitude and the 
mean of the SExtractor uncertainty were calculated for each 
bin. The uncertainty in the colour for each object was deter- 
mined by adding the uncertainty of /Wgis and in quadrature 
(cr^^i = cr^ + cr^). The interval in which 68.3% of the objects 
were closest to this mean colour was used to infer the actual 
Icr colour uncertainty. We assumed that the ratio between the 
old uncertainties cr^ and cr^ was the same for the new un- 
certainties cr'f^ and cr'^. We related between the new and old 
uncertainty in the intermediate and narrowband flux using the 
function crj-j^^^. - '\Ja~ + (^o-fjiter)^, where a is the zero-offset 
for the uncertainty in the flux of bright sources and b is the 
ratio between the new and old uncertainty for the flux of the 
faintest sources. The parameters a and h correspond to imper- 
fections in the photometry and wrongly assumed background 
by SExtractor, respectively. 

Typically, the correction factors are moderate (between 
~30-50%) for the faint sources in the catalogues. Even though 
the correction factors are moderate, we assume that the correc- 
tions for the uncertainties in the broadband B and R are irrele- 
vant, since they are used in a different way than the intermedi- 
ate and narrowband images (see Sect. 14. 2t . 

4.2. Selection criteria 

The following four criteria were applied to select our candidate 
LAEs from the eight initial source catalogues: 
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-0--1 - X " - 



16 '8 20 22 

Fig. 2. Msi5-N colour as seen in the S 1 1 field with the N^w fil- 
ter as the detection image. The mean colour term is ~0.05. The 
heavy bold lines indicate the 68.3% interval of objects colours 
closest to the mean colour in each bin. Each bin contains 200 
data points. The new uncertainty is inferred from this interval. 
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Fig. 4. Galaxy counts as a function of AB-magnitude for the 
A^8i7 filter in the Sll field. Plotted are the Nan source counts of 
the Sll field (crosses) together with the /g^ galaxy counts of 
the Hubble Deep Field (squares). The solid line indicates the 
fitted linear relation between the magnitude and galaxy count. 



1. the narrowband image used as the detection image must 
have the most flux of all the narrowband images and the 
source must have a 4cr detection or better; 

2. the narrowband image with the least flux needs to be a non- 
detection, i.e. less than 2cr; 

3. there must be at least a 2cr detection in the intermediate 
band image; 

4. none of the broadband images, i.e. neither B nor R, must 
have a detection above 2cr. 

Table |2 contains the values of the 2cr detection thresholds of 
the images used for the 6 pixel aperture. In total 33 candi- 
dates were selected using the above criteria. Visual inspection 
showed that 26 sources arose from artefacts of which the vast 
majority were out-of-focus ghost rings from bright stars. The 
final sample contains seven candidate LAEs. 

We note here the importance of the usage of the interme- 
diate band filter If we were to reapply all the criteria except 
for criterion 3, i.e. we do not use the intermediate band images, 
we would obtain 284 candidates instead of the 33 for visual 
inspection. 

The AB-magnitudes, derived line fluxes and luminosities 
for the candidates are shown in Table |6l To convert between 
AZ?-magnitudes and line flux in erg s ' cm"^ we use the follow- 
ing relation: 

Fline = 3 X 10'« 10-0-4('«-+48.590) ^ 

A- 

where AA and A^ are the FWHM and the central wavelength 
of the narrowband filter in A, respectively, and m^g the AB- 
magnitude of the object. In Fig.l^the thumbnails of the seven 
candidate LAEs at z ~ 5.7 are shown. We defer a more detailed 
discussion about the sample properties to Sect.|5l 



4.3. Completeness corrections 

From the Hubble Deep Field (HDF) galaxy number-count data 
for the F814W filter d Williams et al..l996i) we computed com- 
pleteness corrections for our eight source catalogues. The HDF 
counts are determined over the magnitude range /814 = 22 - 29, 
and agree well with our galaxy counts over all narrowband fil- 
ters in the range N -22- 24. FigureEJshows the counts for the 
FS 14IV filter in the HDF and for the Nf,u filter in the S 1 1 field. 
Figurel^also shows the linear fit used as the basis for the calcu- 
lation of the detection completeness. The fit is done to the com- 
bined number count data over two intervals: N^n = [20, 22.5], 
where the WFILAS counts are complete, and /g^ = [22.5,25], 
where the HDF counts are linear. 

Detection completeness is defined as the ratio of WFILAS 
sources to the number expected from the number-count rela- 
tion. Figure |5] shows the derived detection completeness for 
each filter-field combination used for WFILAS. The differ- 
ences are mainly due to unequal exposure times, although fil- 
ter throughput and image quality also play a role. These could 
explain the overall lower sensitivity of the A^824 filter, as can 
be inferred from Fig.|5] Additionally, we correct for detection 
completeness arising due to the intermediate band selection cri- 
terion. We constructed a noise image by stacking the interme- 
diate band images without registering. The completeness is de- 
fined as the rate of recovery of artificially inserted objects. 

Given the different sensitivities of each filter-field combina- 
tion, we define a homogeneous subsample of our initial candi- 
date sample, using the candidates from our four most sensitive 
field-filter combinations. We call this our "complete" sample 
(4 of the 7 LAEs; marked in Table |6j, because once defined, 
we use the curves in Fig. |5]to correct the detected candidate 
numbers for incompleteness, in contrast to our initial "incom- 
plete" sample (all 7 LAEs). The purpose of the subsample is 
that it lies within a uniform flux limit. Figure|5]shows that our 
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Fig. 3. Thumbnails of each region in which the candidate LAEs appears. The thumbnails cover a 19"xl9' 
scale of 0'.'238 pixel"' and North is up and East to the left. From left to right are the filters Mgis, A^gio, N%i-i, 



' region with a 
A^824, B and R. 



E. Westra et al.: WFILAS for Galaxies at z ~ 5.7. II. 



7 




Uncorrected sample 
Completeness corrected sample 



1 


1 1 1 1 1 1 1 1 1 1 1 1 

Nb-d 




^ - " = ^4 

\ 

V \ — 
"\ \ 




\ \ 

\ 

A ^ — 




\ ^ 

V ^ 


: SGP 

1 




22.5 


23.0 ?3.5 7-1.0 

ITIAB - 



Fig. 5. Detection of completeness as function of magnitude de- 
rived from the galaxy density-magnitude relation as described 
in Sect.m From top to bottom the CDFS, SI 1 and SGP fields 
are shown. No N^n data are available for the SGP field. 



four best filter-field combinations consist of the A^gio and A^gn 
filters in both the CDFS and SI 1 fields. These four field-filter 
combinations reach at least 50% completeness at Mab = 23.38, 
or 5.1x10"'^ erg s"' cm"^. We take this as the flux limit of our 
complete sample. As such, the number density derived from the 
complete sample is a more accurate measure of the density of 
sources down to the nominated flux limit than the number den- 



42. S 4:5.0 43.2 43.4 43.6 

Log(Line luminosity (erg s"')) -> 

Fig. 6. Line luminosity distribution of the candidate LAEs. Two 
samples are indicated: all the candidates, but not corrected for 
completeness (solid) and the candidates in the complete sam- 
ple, i.e. candidates of the four deepest narrowband images with 
a magnitude cut-ofF at 50% completeness of the worst of these 
four images (dashed). 

Log(L (ergs"')) Log(cl) (Mpc"' 0.1 dex"')) 



42.7 


-4 83+" -' 


42.8 


-4 65+""* 


42.9 


-4 35+" '3 


43.0 


A 4^+0.14 
^•^-'-0.21 


43.1 


-4 65+" '* 

•"-^-0.30 


43.2 


-5 13+"" 

-^-0.87 


43.3 


c 91+0.19 


43.4 


"-'■■-'^-0.44 



Table 4. Number density of LAEs per luminosity bin as indi- 
cated in Fig. Et. 



sity of the incomplete sample. Figure |6l shows the luminosity 
distribution of the complete sample alongside our initial candi- 
date list, which we call the "incomplete" sample. It shows that 
in using completeness corrections our detected source density 
is up by 50%. 

5. z ~ 5.7 Candidate LAE Catalogue 

In the previous Sect, we introduced two sets of can- 
didate LAEs: the full (but incomplete) sample of seven 
candidate LAEs and a subsample thereof, complete to 
Fum =5. 1x10"'^ ergs"' cm"^ (the complete sample). The flux 
limit of the incomplete sample is almost twice the limit of the 
complete sample (3.4x10"'^ erg s"' cm"^). 

To examine the luminos ity distribution o f our sample we 
use the Schechter function (ISchechtejll976l) . as it is a good 
representation of the data at bright luminosities. From this, the 
luminosity density X of a distribution with a limiting luminos- 
ity Lii„, is given by 



£(L > Ln„,) = fL*Y(a + 2, LuJL*), 



(2) 
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Histograms: 
WFILAS 
Ajikl el al. (2003) 

Function fits: 
WFILAS + Ajiki et al. 
Ajiki et al. only 



42.8 43-0 43.2 43.4 

Log(Line luminosity (erg s"')} -> 

b) ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ 

WFILAS + Ajikl et qI. 

Ajiki et al. only 



the total luminosity density when = 0. We h ave done this 
for three cases: for the candidates of lAiiki et al.l (case A), the 
complete sample of our candidates (case B) and a combined 
sample of these two surveys (case C). For our co mplete sam- 
ple we derive a higher L* (+0.12dex; case B) than lAiiki et alJ 
case A) which implies an increase of the luminosity den- 
sity £. of ~30% . If we scale the luminosity contribution of the 
candidates from lAiitd et al.l to our volume and combine the two 
samples, L* is higher (logL* - 42.66; case C). Table sum- 
marises the results. Detecting LAEs of such bright luminosity 
at this redshift demonstrates the necessity of wide field surveys, 
such as WFILAS, to provide a sample of LAEs at the bright 
end. 

As a second approach, we tried fitting a Sche chter function 
to the combined WFILAS and lAiiki et alJ ll2003l) dataset, using 
a minimised fit (Fig. 17} . We did not use the two lowest lu- 
minosity bins of lAiiki etal.. (.2003i) to constrain the fit because 
these force the function to decline at the faint end. Instead, we 
set the faint end slope to a = - 1.53, similar to the Ha lumi- 
nositv func tion at z ~ 0.24 from lFuiita et al.l (l2003l) . on which 
lAjiki et al.l based their work. FigureEb shows a strong correla- 
tion between L* and (p* due to the slow turn-over at the bright 
end. 



Fig. 7. a) Line luminosity distribution of the complete sample 
of candida te LAEs (solid his tograms) together with the candi- 
dates from lAiiki et al.l ll2003l) (dashed histograms). Both sam- 
ples are corrected for completeness. The errorbars are derived 
using Poisson statistics. Furthermore, two Schechter function 
fits are indicated: one to the combined WFILAS and Ajiki sam- 
ple (dotted) and one to Ajiki sample only (long dashed), b) 
The 68.3%, 95.4% and 99.7% confidence Hmits for the fitting 
parameters L* and (p*. See text for details. 



where a and (p* represent the slope of the faint end of the 
Schechter function and the normalisation constant of the galaxy 
density, respectively. F is the incomplete gamma-function. 
Currently, the luminosity function for LAEs at z ~ 5.7 is poorly 
defined and authors commonly adopt either one or two of the 
three parameters from low redshift surveys to calculate the 
third. 

We examine the influence of non-detections of bright (L > 
L*) LAEs for th e total Lyg luminos ity density by employing the 
same method as lAiiki et al.l ( l20o'3l) . another narrowband imag- 
ing survey aimed at finding LAEs at z ~ 5.7. In the inter- 
est of comparison , we follow lAiiki et akl exactlv and adopt the 
iFuiita et al.lll2003l) values for a (-1.53) and (p* (10^2.62 
Their approach was to solve Eq. (|2} for L*, instead of fitting 
a Schechter function. Fixing <p* and allowing L* and a to vary 
imposes a strong prior on the final fit, it allows us to compare 
directly to the results of lAiiki et al.l bv preserving their method. 
The luminosity density X, was calculated by summing the lu- 
minosity of all candidates (corrected for completeness) and di- 
vided by the corresponding survey volume. With the given sur- 
vey Umits the equation can be solved forL*. Equation (|2} yields 



From the fitting there are three results to conclude. Firstly, 
incorporatin g the four complete ness-corrected WFILAS galax- 
ies into the lAiiki et al galaxies better constrains the 
bright end of the luminosity function. Furthermore, it seems 
that the current generation of surveys is only just reaching the 
volume coverage necessary to discover LAEs with L> L*. The 
histogram in Fig. shows a decreasing number of sources at 
the faint end. At face value, this could suggest that the ionising 
flux of the less luminous sources may be insufficient to escape 
the slowly expanding envelope of neutral hydrogen that sur- 
rounds the H n region in the LAE. Consequently, the sources 
are undetected and the faint end of the luminosity distribution 
decreases. However, it is difficult to detect faint LAEs and so 
the possibility of detection incompleteness cannot be ruled out. 

Figure |8l shows the sky distribution of our candidates in 
each field. All candidates but one are in the CDFS and Sll 
fields. The only candidate in the SGP field is brighter than the 
candidates in the other fields (line flux ~10"'ergs"' cm^^). 
The reason for this is that the M^is filter for the SGP field has a 
shorter exposure time and lower signal-to-noise than the other 
fields. 

In the CDFS field we note that our three candidates ap- 
pear to be spatially clustered. Additi onallv, we not e that the 
confirmed z - 5.78 /-drop galaxy of iBunker et alJ (l2003l) is 
at the same redshift as the WFILAS candidates in this field, 
just like four can didate LAEs from a narrowband survey by 
lAiiki et al.l ( I2OO5I) . We did not detect these four candidates 
since the y are fainter than t he detection limits of WFILAS in 
this field. IWang et ^ il2005ll have also done a narrowband sur- 
vey of the CDFS field. They also find evidence for an overden- 
sity of z ~ 5.7 sources in this field. Similarlv. iMalhotra et alJ 
(I2OO5I) find an overdensity at redshift 5.9 + 0.2 in the HUDF. 
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Fig. 8. Sky distribution of candidate line emitters per field with 
North up and East to the left for the a) CDFS, b) SU and c) 
SGP fields. The "1", "2" and "3" labels correspond to the A^bio, 
A^8i7 and A^824 filters, respectively. The gridlines are separated 
by 7:5. In the CDFS field the Hubble Ultra Deep Field (dashed) 
and GOODS-S (solid) ha ve been indicated. T he confirmed i- 
drop galaxy at z - 5.78 of lBunker et al l(l2003[ cross), LAEs of 
lAiiki et alT ( l2005l plus) and lStanwav et al .I (l2004 triangle) are 
also indicated. In the CDFS field there seems to be a overden- 
si ty of candidates tow ards the southern part of the field, similar 
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Fig. 9. A 38" X 38" region around the confirmed LAE in the 
SGP field. The image is created from the pre-image taken with 
VLT/FORS2. The image has a pixel scale of 0'.'252pix North 
is up and East is to the left. 

6. Confirmed LAEs^ 

In IWestra et alJ (l2005l) we reported the spectroscopic follow- 
up of one of the candidates, Jl 14334.98-014433.7 (Sll_13368 
in that paper, hereafter Sll_5236''). It was confirmed to be a 
LAE at z = 5.721. Here we present the spectral confirmation of 
a new candidate, J004525.38-292402.8 (hereafter SGP.8884), 
at z = 5.652. We also show its pre-imaging and compare its Lya 
profile to S 11.5236. SGP_8884 and SI 1.5236 are the only two 
out of the seven candidates presented in this paper for which 
we have obtained spectra. 

6.1. Spectral data reduction 

A pre-image with an intermediate band filter (FWHM =13 nm) 
centred at 8 15 nm was taken with VLT/FORS2 on 2005 August 
9. The 0'.'252pix ' plate scale undersamples the <0.5"stellar 
point spread function of the frames which were taken during 
excellent seeing. SGP_8884 is unresolved, implying that the 
FWHM of the emitting region is <2.2 kpc. A 38" x 38" region 
around the object is shown in Fig.|5] 

The spectroscopy consists of four exposures of 900 s, taken 
on 2005 October 3 with FORS2 using the 1028z grism and a 
1" slit. The frames were overscan subtracted and flatfielded. 
They were combined by summing individual frames, thereby 
removing cosmic rays in the process. 



' Based on observations made with ESO Telescopes at the Paranal 
Observatories under programmes ID 076.A-0553 and 272.A-5029. 

* The object names are derive d from SExtractor I Ds. Refinements 
to our detection procedures since IWestraetalJ<2005t) caused a change 
in the ID and, therefore, in the object name 



10 



E. Westra et al.: WFILAS for Galaxies at z ~ 5.7. II. 



Velocity {km s ') ^ 
-1000 1000 2000 
1.5x10-"! 




8060 8080 8100 8120 81^0 

Wavelength (A) ^ 



Fig. 10. (Top) Flux calibrated spectrum of the confirmed can- 
didate LAE SGP_8884, the brightest candidate in our sample. 
The histogram shows the observed spectrum. Indicated in grey 
is the best-fitting single component model after convolution 
with the instrumental profile. The heavy bold lines indicate 
three regions for which we have calculated a mean continuum. 
(Bottom) Residuals from the observed data minus model fit. 
The histograms indicate the Icr-error spectrum from the ob- 
served data, which includes both sky- and Poisson noise. The 
feature at 8125 A is due to a remnant cosmic ray from one of 
the spectral frames. 



The spectrum was flux calibrated using a standard star 
(HD 49798) taken with a 5" slit and corrected for slit-loss. 
This was calculated assuming a Gaussian source profile with 
a FWHM of 0'.'72 as measured from the spatial direction of the 
spectrum. The flux lost due to the 1" slit was calculated and 
added to the spectrum of the object. 

6.2. Line fitting 

Figure [TO] shows the reduced spectrum of SGP_8884 along- 
side its best model fit. The spectrum has an asymmetric line 

S irofile, similar to our previously confirmed candidate LAE 
Westra et al.ll2005l) . It unlikely originates from a redshifted 
[On] line at z ~ 1.2 because the resolution of our spectrum 
is high enough to resolve the [On] /1/13726,3728. Figure fTTI 
shows the spectrum of one such [On] emitter at z = 1.18 
which was included in the same observations as SGP_8884. 
Furthermore, we do not find any other spectral features in our 
spectrum, such as Hfi or [Nn], which could classify the emis- 
sion coming from a lower redshift galaxy. Hence, we identify 
the line as Lyo' at z = 5.652. With a total spectral line flux of 
(1.0 + 0.1) X 10"''' ergs"' cm"^ (slit-loss corrected), SGP_8884 
is the brightest LAE at redshift -5.7 to date. The line flux 
derived from the spectrum is consistent with the flux derived 
from narrowband photometry (9.5 + 1.4) x 10"'^ erg s"' cm"^, 
which is given in Table|6] The spectral line flux corresponds to 
a line luminosity of Liine = 3.5 x lO'*-' erg s"' and a star forma- 
ti on rate of 32 Mq yr " ' , using the star formation conversion rate 
of lAiiki et al.1 (l2003h . If we adopt ~16 pixels (= 32kpc^) as an 
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Fig. 11. (Left) Two dimensional spectrum of an [On] emitting 
galaxy at z = 1.18 in the SGP field. (Right) The extracted 
one dimensional spectrum. We are able to easily resolve the 
[On] /l/i3726,3728 lines with the available resolution. Both 
spectra are background subtracted. 

upper limit to the size of the emitting region, we derive a star 
formation rate surface density of > 1 Mp yr"' kpc" ^. 

Following ea rlier works (e.g. lPawson et a l.*2002 UHu et alJ 
l2004tlWestra et al.ii2005i) we fitted a single component model 
to the Lya line SGP_8884. The model consists of a truncated 
Gaussian with complete absorption blueward of the Lycf line 
centre. We find an excess of flux in the observed data com- 
pared to the model around 81 10 A. This suggests the pres- 
ence of a second line component redward of the main peak. 
To test this, we measured the mean continuum levels, both 
red- and blueward of the line, as well as across the red- 
flanking region of the line. The continuum is calculated as 
the weighted mean of the flux density over this region. This 
yields for continuum in the red-flanking region a flux density of 
(3.2 + 0.8) X lO"'"* erg s"' cm^^ A"' . Red- and blueward of the 
Lycf line the continuum is (-1.0 + 0.8) x 10"'^ erg s"' cm"^ A"' 
and (0.9 ± 0.6) X 10"'^ erg s"' cm"^ A"', respectively. These 
continuum levels are indicated by the heavy bold lines in 
Fig. [TO] The lower limit for the rest frame equivalent width 
derived from the continuum of the red flank is 46 A. The rest 
frame equivalent width derived from the 2cr upper limit of the 
continuum redward of the line is 125 A. 

To see if the excess of flux in the red flank of the 
Lya line can be explained by an outflow, we fit a sec- 
ond Gaussian component to the spectr um of SGP_8884, a s 
we did to the spectrum of Sll_5236 in IWestra et alJ ( l2005l) . 
This yields an extremely faint and broad second component 
(/peak ~5x 10"''' erg s"! cm"- A"' and FW//M~1700kms"'). 
The precise parameters for the red component are difficult to 
constrain given its faint and broad profile. The parameters from 
the single component model for SGP_8884 and the single and 
double component models for SI 1_5236 are given in Table|5] 

6.3. Discussion/Comparison 

The Lya emission we see is due to intense star formation 
rates synonymous with local starburst galaxies. Star forma- 
tion rates per unit area in excess of 0.1 Mq yr"' kpc^ are prone 
to produce large scale outflows of neutral hydrogen from a 
galaxy, po wered by the sup ernovae and stellar winds of mas- 
sive stars ( iHeckmanlEooa . The most efficient way for Lyor 
to escape from the compact star forming regions is due to 
scattering of the photo ns by the entrained neutral hydrogen 
dChen & Neufeldll994 . The kinematics and orientation of the 
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Fig. 12. Comparison of the Lya line profiles of the two WFILAS sources, SGP_8884 and Sll_5236. a) Two dimensional 
background-subtracted spectrum of SGP_8884. b) Observed Lya line (histograms) with the best-fitting one component model 
(grey solid line), c) Same model line profile as in b) but before convolution with the instrument profile, d) Observed data minus 
model fit (as plotted in b)). Also shown (histograms) is the Icr-error spectrum from the observed data, which includes both sky- 
and Poisson noise. Panels e) through h) show the same for Sll_5236. The horizontal axes show both the wavelength (in A; 
bottom) and velocity offset from the centre of the full Gaussian of the Lya line (in km s ' ; top) 



outflowing neutral hydrogen can alter the Lya profile by ab- 
sorbing photons bluer if along the line of sight, or backscatter- 
ing re dder than Lyo' if behind and receding (e.g. lPawson et al] 
I2OO2I) . Lya emission can also arise when large scale shocks 
from starburst winds impinge on clumps (~100pc) of con- 
dense d gas accreting onto the halo jBland-Hawthorn & Nulser] 
l2004 . 

Most examples of asymmetric Lya emission at z ~ 6 show 
an extended tail implying backscattering over a fairly wide 
ra nge of velocities beyond the central Lya emission (e.g. Fig. 9 
of lHuetalJf2004l) . The Hmiting physical size of SGP_8884 
(FWHM <2.2kpc) is consistent with the scale of emitting re- 
gions in the local starburst galaxy M82 which span 0.5 to 
Ikpc dCourvoisier et alJll99(i iBlecha et alJll990l) . This, and 
the scale of its outflow, make it fairly typical of both the star- 
bursting sources seen at z ~ 6 and their local counterparts. 

The tentative dis covery of a second component in S 11 _5236 
dWestra et al]l2005h could be explained by either an expand- 
ing shell of neutral hydrogen jPawson et alJl2002l lAhn et all 



'20031, or by infall of the IGM onto the LAB jPiikstra et alJ 
120051) . The flux of the intrinsic Lyo' line depends heavily on the 
model. It is suggested that the total intrinsic Lya flux emerging 
fro m these sources is underestimated by an order of magnitude 
(e.g. lPiikstra et alJ) . Therefore, the star formation rates derived 
from the observed Lya lines could be heavily underestimated. 

Figure [12] shows a comparison between the line profiles 
of the two LAEs discovered with WFILAS. Sll_5236 differs 
from SGP_8884 in that a clear pe ak, -20- 90km s~ ' wide, is 
seen -400 km s'^ redward of Lya dWestra et alJ200-'^ . The red 
component is narrower (~15%) and relatively stronger than 
SGP_8884. The difference in the width of the red component 
is even more pronounced (~30%) when we compare the main 
peak of the two-component fits to the spectrum of Sll_5236 
to the single peak of the one-component fits to the spectrum of 
SGP_8884. This can clearly be seen in panels a and e of Fig.fT^ 

Ultimately, such outflows are thought to be responsi- 
ble for the chemic al enrichment of the IGM by z ~ 6 
dAguirre et alJ200ll) . Outflows are a process facilitating the es- 
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Component /peak FWHM Av 



(1) (2) (3) (4) (5) (6) (7) 



SGP_8884 single component 

Single peak 8086.2 1.2 xlO"" 15.7 


580 






SI 1_5236 single component 

Single peak 8172.2 8.3x 10-"* 13.5 


495 






S 11-5236 double component, "broad" 
Main peak 8173.1 8.0 xlO""* 11.3 
Red peak 8184.1 1.9 xlO^"* 2.3 


413 
85 


+ 11.1 


+406 


S 11-5236 double component, "narrow" 
Main peak 8173.1 8.1 xlO""* 11.2 
Red peak 8184.1 4.8 xlO^"* 0.5 


413 
18 


+ 11.0 


+403 



Table 5. Parameters for the single component model to 
SGP_8884 before convolution with the instrumental profile. 
We also include the parameters for the single component and 
the two double c omponent models o f the previously confirmed 
LAE Sll_52 36 j Westra et all f2005l) . These parameters differ 
slightly from IWestra et all since we have subsequently cor- 
rected the spectrum of S 11 -5236 for slit-losses. 
Notes: (1) component of the fit, (2) central wavelength 
of the fitted component in A, (3) peak flux density in 
ergs"' cm A-\ (4) and (5) FWHM of fuU Gaussian of the 
profile in A and kms"', respectively, (6) and (7) velocity shift 
of the second component in A and km s"' . 



cape of UV photons, which are the origin for the UV back- 
ground (fMadau et al..l99^) . 

7. Summary 

In this paper we have presented the Wide Field Imager 
Lyman- Alpha Search (WFILAS), which uses a combination of 
narrow-, intermediate and broadband filters on the ESO/MPI 
2.2m telescope to search for LAEs at redshift z ~ 5.7. This 
search has resulted in seven bright (L > 1 . 1 x 10"*'' erg s"' ) candi- 
date galaxies across three fields spanning almost 0.8 sq. degree. 

Most of our candidates are in the regimes of bright lumi- 
nosities, beyond the reach of less voluminous surveys. Adding 
our candidates to those of earlier such surveys results in an 
integrated luminosity density X ~30% higher than found by 
such surveys alone. We also find potential clustering in our 
CDFS field, supporting overdensities discovered by other sur- 
veys. Spectroscopic follow-up for confirmation in this area will 
be crucial. 

Two candidates have been confirmed to be LAEs at z ~ 5.7 
by means of spectroscopy. One of these galaxies is the brightest 
LAEs at this redshift. The broad, asymmetric profiles of the 
Lya line of both objects are consistent with neutral hydrogen 
backscattering of a central starbursting source. 
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SExtractor ID 


Object ID 


B 


R 




A'sio 


A'817 


A'824 


Line flux 


Luminosity 




















(10""ergs-'cm-2) 


(10-*^^ ergs-') 


CDFS_1864'' 


J033215. 14-280013.9 


>26.25 


>26.56 


24.72 ± 


0.46 


23.14 ± 0.26 


>24.27 


>23.93 


6.5 + 1.5 


2.3 ± 0.5 


CDFS_4928" 


J033 145. 97-2753 16.4 


>26.25 


>26.56 


24.59 ± 


0.41 


23.38 ± 0.32 


24.11 +0.47 


23.61 ± 0.41 


5.2 + 1.6 


1.8 ±0.5 


CDFS-SSSS" 


J033202.37-275211.3 


>26.25 


>26.56 


24.70 ± 


0.45 


23.32 ± 0.31 


>24.27 


>23.93 


5.5 + 1.5 


1.9 ±0.5 


SI 1.5236"'' 


Jl 14334.98-014433.7 


>26.63 


>26.59 


24.31 ± 


0.42 


>24.13 


23.05 ±0.18 


>23.74 


7.0 + 1.2 


2.5 ± 0.4 


S11.892r 


Jl 14218.90-013544.6 


>26.63 


26.38 ± 0.45 


23.88 ± 


0.28 


23.98 ± 0.47 


23.41 ± 0.26 


>23.75 


5.0 + 1.2 


1.8 ± 0.4 


Sll_10595 


Jl 14312.46-013049.6 


>26.63 


>26.60 


24.44 ± 


0.47 


>24.13 


23.52 ± 0.28 


>23.75 


4.5 ± 1.2 


1.6 ± 0.4 


SGP.8884'' 


J004525.38-292402.8 


>26.07 


>26.41 


23.33 ± 


0.20 


22.73 ± 0.16 




>24.06 


9.5 + 1.4 


3.3 ± 0.5 



Table 6. The candidate list of the WFILAS survey after the selection as described in Sect.|4] From left to right are the object 
name, the B, R, M^is, A^gio, A^gi? and A^s24 AB-magnitudes, line flux calculated from the narrowband magnitude in which the 
object was detected and line luminosity. For all measurements less than 2cr the 2cr upper limit has been given. 



" Galaxy is in the complete sample 

* Confirmed LAEatz = 5.721. See text and lWestra et all llOOj) for details. 
Signal-to-noise in the range 2 - 3(T for R band in the 10 pixel aperture, but < 2cr in the 6 pixel aperture 
Confirmed LAE at z = 5.652. See text for details. 





a 


log0* 
Mpc-^^ 


logL* 
erg s"' 


log L;,„, 

erg s"' 


log V 
Mpc' 


log£ 

erg s"' Mpc"'' 


Comment 


Case A 








42.85 


5.26 


39.04 


Sum of the candidates from Aiiki et al. (2003) 




-1.53 


-2.62 


42.61 


42.85 


5.26 


39.04 


Inteerated luminositv function down to Aiiki et al. (2003) survey limit (7.0x10"'^ ere s"') 




-1.53 


-2.62 


42.61 






40.27 


Integration of the entire luminosity function 


CaseB 








43.26 


5.71 


38.36 


Sum of the candidates from completeness corrected WFILAS sample 




-1.53 


-2.62 


42.74 


43.26 


5.71 


38.36 


Integrated luminosity function down to the limit of the completeness corrected sample (1.8x10'" erg s"') 




-1.53 


-2.62 


42.74 






40.39 


Integration of the entire luminosity function 


CaseC 








42.85 


5.84 


39.19 


Sum of the combined WFILAS and Aiiki et al. (2003) samples low luminositv corrections 




-1.53 


-2.62 


42.66 


42.85 


5.84 


39.19 


Intesrated luminositv function down to Aiiki et al. (2003) survev limit (7.0x10''^ era s"') 




-1.53 


-2.62 


42.66 






40.32 


Integration of the entire luminosity function 



Table 7. Calculation of the Schechter function parameter L* and luminosity density X according to Aiiki et al. ( 2002J for their sample, our complete sample and the combination 
of the two. For each sample the luminos itv density has bee n derived from the sum of the candidate luminosities divided by the corresponding survey volume. Then Eq. Q was 
solved for L*, with given a and cf>' from lAiikiet"al]ll2003l) . Finally, the entire luminosity function was integrated to give the final luminosity density. 
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List of Objects 



Jl 14334.98-014433.7 on page 9 



J004525.38-292402.8t on page 9 



